Nernst effect and diamagnetism in phase fluctuating superconductors 
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When a superconductor is warmed above its critical temperature T c , long range order is de- 
stroyed by fluctuations in the order parameter. These fluctuations can be probed by measurements 
of conduct ivity[l[, diamagnetism 0] and of the Nernst effect 0, 0, [1| . Here, we study a regime where 
superconductivity is destroyed by phase fluctuations arising from a dilute liquid of mobile vortices. 
We find that the Nernst effect and diamagnetic response differ significantly from Gaussian fluctu- 
ations - in particular, a much sharper decay with temperature is obtained. We predict a rapid 
onset of Nernst signal at a temperature T onse t that tracks T c , rather than the pairing temperature. 
We also predict a close quantitative connection with diamagnetism - the ratio of magnetization 
to transverse thermoelectric conductivity a xy reaches a universal value at high temperatures. We 
interpret Nernst effect measurements on the underdoped cuprates in terms of a dilute vortex liquid 
over a wide temperature range above T c . 



In recent years, the Nernst effect has emerged as an 
important probe of strongly-correlated electron systems. 
The Nernst signal is the electric field (E y ) response to a 
transverse temperature gradient V X T ', 
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A large Nernst signal has been detected in the nor- 
mal state of quasi-2d cuprate[|| and heavy fermionQ 
samples, and in thin filmsQ, well above T c . This is 
in contrast with typical (non-ambipolar) metals.where 
the Nernst effect is usually weak: it was shown[6| that 
Fermi liquid quasiparticles with energy-independent scat- 
tering rates do not contribute to the Nernst signal. This, 
together with the proximity of the large-Nernst region 
in some of these materials to superconducting phases, 
points towards fluctuating superconductivity as one nat- 
ural source for the Nernst signal. 

Theoretical studies of the Nernst effect in cuprate su- 
perconductors include the analysis of Gaussian fluctu- 
ations above the mean-field transition temperature (Io|. 
and a Ginsburg-Landau model with interactions between 
fluctuations of the order parameter In this Letter, we 
consider the Nernst effect due to thermal fluctuations in 
a phase-only model. Throughout, we assume that the 
superconducting order parameter ip(x) = Aoe t9 ^ has a 
frozen amplitude Ao(x) = const. Such a situation can 
arise in granular thin films and in Josephson junction ar- 
rays, where the order parameter on individual grains is 
well-established, whereas T c is given by the weak Joseph- 
son coupling between grains. This picture may also be 
appropriate for the underdoped cuprates, where the pair- 
ing gap is thought to be much larger than the transition 
temperature, ksT c <C A0, Q. Here, superconductivity 
is destroyed by loss of phase coherence via thermally gen- 
erated vortex anti-vortex pairs. Phase fluctuations from 



vortex diffusion have been proposed as the dominant con- 
tribution to the Nernst signal 3, 9]. This vortex picture is 
most useful in the dilute limit, when the spacing between 
vortices (both field induced and thermally generated) is 
much larger than their core radius (the zero temperature 
coherence length £q)- Then, vortices have a well defined 
identity and the amplitude is suppressed only in the small 
area of the sample occupied by vortex cores. 

Our starting point is the Lawrence-Doniach model of 
a layered superconductor, 
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where i,j label lattice points within a layer, and n,m 
label the layers. The lattice vector potential due to an 
external magnetic field, Aij = =£■ J r 3 dr ■ A, is static and 
unscreened, corresponding to an extreme type-II super- 
conductor. We consider the limit deep in the ordered 
phase within mean field theory — r 3> ksT, where phase 
fluctuations dominate r/>,- )7l = Aoe* i,n . This reduces the 
model above to an XY model. The inter-layer coupling 
J±_ stabilizes true long-range superconductivity. How- 
ever, we have verified that realistic values of J71 increase 
T c relative to the 2D Kosterlitz-Thouless transition Tkt 
only by a small amount, and that the normal state prop- 
erties of interest are not significantly affected by J\_ , ex- 
cept very close to T c . Hence, in what follows we set 
J±_ = and consider the 2D XY model with Josephson 
coupling J = A 2 J. 

To study transport, we supplement this statistical me- 
chanics model with model-A Langevin dynamics, corre- 
sponding to interaction with a heat bath that does not 
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preserve any conservation laws 14| : 
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Here, r provides a characteristic time scale for the dy- 
namics. The stochastic noise rji{t) is Gaussian corre- 
lated, with a variance chosen to satisfy the fluctuation- 
dissipation theorem, 



(Vi(.t)vj(t')) = 2k B T T 5 lJ 5{t-t'). 
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The model ([3]), and (0J has only three free parameters: 
J, r, and the lattice constant, a (or equivalently a field 
scale, i/o = < I ) o/(27ra 2 ), defined via the superconduct- 
ing flux quantum <& ). Of these, J is an overall energy 
scale, set by fixing Tkt- The length scale a depends on 
the physical system in question. For Josephson junction 
arrays or granular superconductors, this is the spacing 
between grains. For uniform superconductors, a can be 
determined by comparing the correlation length in the 
XY model away from the transition, eg. at T = 2T C , 
with the typical separation between thermally induced 
vortices at that temperature. Thus, a is given by a com- 
bination of vortex fugacity and core radius £ [13] ■ Within 
£o, the superconducting amplitude is significantly sup- 
pressed. Hence, the dilute limit, where the separation 
between vortices exceeds £o , determines the temperature 
window over which a phase-only description is appropri- 
ate. Finally, the time scale r does not affect thermody- 
namic quantities, such as magnetization, nor does it enter 
the transverse thermoelectric conductivity, a xy , which is 
closely related to the Nernst effect, 
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where o~ xx is the electrical conductivity, and we have as- 
sumed particle-hole symmetry. Hence, a xy and magne- 
tization predicted by this model are only functions of 
T/T KT and H/H Q . 

By an Onsager relation, the transverse thermoelectric 
conductivity a xy can be obtained either from the elec- 
tric current response to a temperature gradient, jt. r , x — 
—a xy V y T, or from the heat current response to an elec- 



tric field, 
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A third method to compute 



a xy is through a Kubo formula involving the unequal 
time correlator (jx(t)jy(0))- As a check of our numerics, 
and of the use of proper magnetization subtractions 13[ , 
we confirmed that all three methods agree for a repre- 
sentative set of temperatures and magnetic fields. 

Intermediate Temperature Regime: The results for o? x d v 
on a single layer are shown in Figs. [1] and O As noted 
above, a xy is independent of the parameter r. For a 
single layer, its value only depends on the ratios T /Tkt 
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FIG. 1: Transverse thermoelectric conductivity for a single 
plane, in units of the quantum of thermoelectric conductance 
2eksjh. Simulations on a cylindrical geometry, with system 
size ranging from 60 x 60 to 200 x 200. Inset: Diamagnetic 
response. Temperatures as in the main figure. 



and H/Hq, and is expressed naturally in units of the 
2D "quantum of thermoelectric conductance", 2ek B /h. 
To compare with thin film and layered systems, one 
must divide by the film thickness/layer separation d, 
Oi X y = c^y/d- The inset of Fig. [T] shows the diamagnetic 
response. Below T c , the magnetization diverges logarith- 
mically in H. The magnetization edge currents corre- 
spond to a "depletion layer" near the edge of the sample, 
where vortex density is smaller than the density in the 
bulk, H / $o- In the vortex picture, magnetization is anal- 
ogous to the work function in a metal, and we find 17| , 
to leading order in H, 
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This is similar to a previous result [l5(, and is in good 
agreement with our simulations. 

High Temperature Expansion: For T« J, the phase- 
only model allows for an analytically tractable regime 
that is entirely different from the Gaussian regime con- 
sidered previously fioj]. The high temperature expansion, 
carried out in powers of J/ksT, is conveniently per- 
formed using the Martin-Siggia-Rose formalism la. 17 1. 



Since both M z and a xy require a magnetic field, the ex- 
pansion of these quantities involves graphs enclosing fi- 
nite a magnetic flux. The leading term thus depends 
on the smallest closed graph - on a square lattice this 
involves 4 links, and is hence proportional to (J/T) 4 , 
whereas on a triangular lattice it goes as (J/T) 3 : 



a 



M 



2d 

xy 

2d 



X- 



2ek 



= -2A 



J 

h V T 

2ek B (J 
T 



" . H 
sin ■ 



#o 
" . H 

silj 



(7) 



H 



3 



X 
X 




FIG. 2: Contour plot of log [a 2 x d y / {2ek B /h)]. For small T 
and H, an inter-layer Josephson coupling J±, absent in these 
simulations, stabilizes 3d superconductivity. Note that, unlike 
the Nernst signal ejv = a xy /a xx in Ref. f3|, there is no ridge 
field in our simulations of a xy . The difference is likely due to 
the diverging electrical conductivity a xx as T c is approached. 
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Here, // = 4 and A = 7r/8(/x = 3 and A = tt/4) for a 
square (triangular) lattice. Despite the lattice-dependent 
behavior of a xy and M z /T, their ratio ([5]) is equal to 
-2, independent of the lattice. The same value of - 
2 is obtained in the Gaussian regime of the dynamical 
Ginzburg-Landau equationfiol. Hlj . and seems to be a 
robust feature of fluctuating superconductivity at high 
temperatures. The ratio \M Z \/Ta xy , shown in Fig. [31 
is only weakly field-dependent, and tends to 2 at high- 
temperatures. This points to a close quantitative con- 
nection between the Nernst effect and diamagnetism Q ■ 
Comparison to Nernst Measurements in the Cuprates: 
Experimental measurements of both Nernst voltage and 
conductivity are required to obtain a xy . Such experi- 
mental data is available on underdoped La2-sSr x Cu04 
(x = 0.12 and T c = 28 K) in weak fields. This is shown in 
Fig. 21 which displays the Nernst coefficient times con- 
ductivity, va xx — ^ttt ■ To compare, we choose 
in our simulation, J — Ji = 30.2 K, corresponding to 
T c = IMTrt, and H = 50 T (for this sample, H c2 wlOO 
T, hence H C 2 > Hq, consistent with a relatively dilute 
vortex liquid). With these values, the simulation gives 
good agreement with absolute experimental values in the 
regime T c < T < 2T C K, except for the very lowest tem- 
perature point T — 30 K. This is very close to T c , so that 
3d superconducting fluctuations, ignored here, are likely 
dominant. 

Comparison with High Temperature Data: The inset 
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FIG. 3: The dimensionless ratio \M Z \/Ta xy , as a function 
of magnetic field for various temperatures. In the high- 
temperature limit, this dimensionless ratio is expected to sat- 
urate at a value of 2 for all magnetic fields. Inset: \M Z \/Ta xy 
for H = 0.3Lffo vs. T/T KT 



of Fig. Q] shows the measured vo xx on a log-log plot ex- 
tending to T = 1207^ 4T C . The data displays a rapid 
decay over a large temperature range, in general agree- 
ment with our expectations. In particular, a T~ 4 decay 
is observed, which is the high temperature result ((7|) on 
a square lattice. However, since in the high tempera- 
ture regime the precise power depends strongly on lat- 
tice geometry (in contrast to the intermediate temper- 
ature regime) justification for using a nearest neighbor 
square lattice model (as opposed to, say, a triangular 
lattice model) is required. While the underlying d-wave 
symmetry of the cuprates, coupled with the fact that a is 
a microscopic length, about 6 times the lattice spacing, 
may be invoked, whether these are sufficient to justify the 
square lattice model is unclear at present and requires 
further work. 

The characteristic scale of a xy at these temperatures 
is enhanced from what one would expect for a super- 
conductor with T c — 28K. For example, the best fit 
to the square lattice high temperature expansion (solid 
line) requires J = Jh = 52K, larger than the effective 
coupling J i = 30.2 K that yiels the correct T c . This 
may be naturally attributed to thermal d-wave quasi- 
particles, omitted in this analysis, which suppress the 
long distance superfluid density [19( but not the super- 
fluid density at shorter scales [17} , which controls the high 
temperature behavior. The ratio Jh/ Ji = 1.6 is consis- 
tent with measurements of the temperature dependent 
superfluid density in other cuprates (l8| . A prediction 
from this scenario is that magnetization should continue 
to track a xy T. 

Onset Temperature: Ong and collaborators Q define a 
temperature r onsot where the fluctuating contribution to 



4 



H 

E 
01 

> 




0.5 



0.0001 



30 40 60 80 100120 
T(K) 



30 



35 40 
T(K) 



45 



50 



FIG. 4: Comparison to underdoped La2-a;Sr T Cu04 (x = 
0.12, T c = 28 K). The experimental data@|, shifted by 
a constant quasiparticle background contribution vb<J X x = 
— 0.011VV KQTm, is indicated by •. Simulation results are 
shown by X with error bars. The value of the experimen- 
tal point at T = 30K (■) has been divided by 4 to fit in 
the figure. This point is very close to T c , and hence va xx 
is likely dominated by 3d fluctuations, not considered in our 
simulations. 



As a benchmark, we note that BCS theory predicts the 
value tbcs ~ 0.7h. This yields a fluctuating conductivity 
at T = 2T C that is only about 10% of the quasiparticle 
conductivity in this material. 

In conclusion, we have studied the transverse thermo- 
electric conductivity a xy and the diamagnetic response 
M z in the classical XY model with model-A dynam- 
ics. We have obtained numerical results at low tempera- 
tures, and analytic results at high temperatures, that are 
functions only of two variables, T/T c and H/Hq, where 
H < H c2 is a characteristic field scale set by vortex pa- 
rameters. In our model, we predict that a xy and M z for 
different systems (e.g. different dopings) should collapse 
into a single curve when expressed in terms of the system- 
dependent T c and Ho. We show that M z /T and a xy track 
each other and, in particular, we predict that their ratio 
tends to —2 at high temperatures. Measurements of a xy 
on the underdoped cuprate Lai.8sSro.i2Cu04 display a 
sharp temperature decay, in agreement with our model. 

We would like to thank PW. Anderson, D. Huse, S. 
Kivelson, J.E. Moore, S. Mukerjee, N.P. Ong, and Y. 
Wang for many insightful discussions. 



the Nernst effect can no longer be experimentally distin- 
guished from the quasiparticle background. Here, since 
we have a natural scale for a xy , we define T onsot as the 
temperature where a xy has decayed to a small fraction 5 
of the quantum of thermoelectric conductivity, 

2ek H . 
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For our model, a 2 x d y = ^-F (T/T KT , H/H ), hence 
Tonsot is proportional to T c . The essential point is that, 
because a xy depends strongly on temperature, when in- 
verted, T onset is only a weak function of 6 and H. For 
instance, the choice H = Hq /4 and 6 = 0.01 on the 
square lattice yields 
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This is consistent with the observation that the 
experimentally-defined T onset roughly tracks T c as doping 
is varied . Due to the strong temperature dependence 
of a xy , the onset of Nernst effect is very sharp, in con- 
trast to Gaussian fluctuations [Io| , where a xy only decays 
as 1/(T C — T) at high temperatures. 

Experimentally, measurements of the electrical con- 
ductivity a xx do not have discernable contributions due 
to fluctuating superconductivity at temperatures of or- 
der T onset [H . Since a xx is proportional to the parameter 
r appearing in Eq. Q, this places a constraint on the 
maximum value of r. The high temperature expansion 
yields 
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